function S=FSimpson(f,a,~,N)
%£表示被积函数句柄
% a，b表示被积区间[a，b]的端点
%N表示区间个数
%S是用复化Simpson公式求得的积分值h=(b-a)/N;
fa=feval(f,a);
fb=feval(f,b);
S=fb+fa;
x=a;
for i=1:N
x=x+h/2;
[x=feval(f.x);
S-S+4fx;
x=x+h/2;
[x=feval(f,x);
S=S+2*[x;
end
S=hS/6;